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ABSTRACT 

VLB A high spatial resolution observations of the disk structure of the active 
galactic nucleus NGC 1068 has recently revealed that the kinematics and geom- 
etry of this AGN is well characterized by an outer disk of H2O maser emission 
having a compact milliarcsecond (parsec) scale structure, which is encircling a 
thin rotating inner disk surrounding a ~10'^ Mq compact mass, likely a black 
hole. A curious feature in this source is the occurrence of a misalignment be- 
tween the inner and outer parts of the disk, with the galaxy's radio jet being 
orthogonal to the inner disk. We interpret this peculiar configuration as due to 
the Bardeen-Petterson effect, a general relativistic effect that warps an initially 
inclined (to the black hole equator) viscous disk, and drives the angular momen- 
tum vector of its inner part into alignment with the rotating black hole spin. We 
estimate the time-scale for both angular momenta to get aligned as a function the 
spin parameter of the Kerr black hole. We also reproduce the shape of the parsec 
and kiloparsec scale jets, assuming a model in which the jet is precessing with 
a period and aperture angle that decrease exponentially with time, as expected 
from the Bardeen-Petterson effect. 
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1. Introduction 

Recent mid-infrared interferometric observations (Jaffe et al. 2004) were able to resolve 
the inner dust torus of the nearby (14.4 Mpc; Bland- Hawthorn et al. 1997) Seyfert 2 galaxy 
NGC 1068. The emission distribution along this torus showed a small hot source (probably 
the accretion disk) , surrounded by a parsec-scale warm dust structure. 

The hot source coincides with the bright H2O masers (Greenhill et al. 1996; GaUimore et 
al. 1996b, 2001, 2004) detected with Very Long Basehne Interferometry (VLBI). The maser 
emission seems to originate in a sub-Keplerian rotating disk with inner and outer radius 
corresponding to ~0.65 and 1.1 pc respectively enclosing, within its inner radius, a black 
hole and an accretion disk of ~ 10^ Mq. The departure from Keplerian motion has been 
attributed to the high mass of the accretion disk that dilutes the gravitational field of the 
black hole (e.g., Hure 2002; Lodato & Bertin 2003). 

The continuum emission at sub-arcsecond scales revealed four components arranged in a 
jet-countcrjet structure, with the core coinciding with the accretion disk position (Wilson & 
Ulvestad 1987; Muxlow et al. 1996; Gallimore et al. 1996a, 2004). The radio jet orientation 
changes monotonically with its distance from the core, bending towards the northeast direc- 
tion. At larger scales, extended lobes are observed, with the northeast lobe having a conic 
shape, similar to the bow shock structure produced by jet-environment interaction (Wilson 
& Ulvestad 1987). 

An interesting result from maser observations is that the maser disk is not completely 
aligned with the major axis of the accretion disk across all its extension (e.g., Greenhill 
& Gwinn 1997; GaUimore et al. 2004). Indeed, they reveal a maser distribution similar to 
a warped disk, with the inner parts almost perpendicular to the jet, but deviating from 
this configuration as the distance to the core increases. Gallimore et al. (2004) proposed a 
configuration in which the thin hot disk is perpendicular to the radio jet, while the misaligned 
maser disk points towards the parsec-scale dusty-torus. 

In this work we investigate the possibility that both the bent jet and warped accretion 
disk are a consequence of the Bardeen-Petterson effect, which results from the misalignment 
between the angular momenta of the Kerr black hole and disk. In § 2, we discuss the 
Bardeen-Petterson mechanism, used to calculate the time-scale for alignment between the 
two angular momenta. In § 3 we present the physical parameters of the accretion disk of 
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NGC 1068, as constrained by the observations found in the hterature, and we use them to 
model the warping of the accretion disk due to the Bardeen-Petterson effect as a function of 
the ahgnment time-scale. In § 4, we show how jet morphology can be used independently to 
estimate the alignment time-scale. Conclusions are presented in § 5. 



Frame dragging produced by a Kerr black hole causes precession of a particle if its 
orbital plane is inclined in relation to the equatorial plane of the black hole. The precession 
angular velocity Olt due to the Lense-Thirring effect is given by (e.g., Wilkins 1972): 



where G is the gravitational constant, c is the light speed, r is the radial distance from a 
black hole of mass Mbh and angular momentum Jbh, defined as Jbh = a*GM|jj/c, with a* 
the ratio between the actual angular momentum and its maximum possible value. 

The combined action of the Lense-Thirring effect and the internal viscosity of the accre- 
tion disk forces the alignment between the angular momenta of the Kerr black hole and the 
accretion disk. This is known as the Bardeen-Petterson effect (Bardeen & Petterson 1975) 
and affects only the innermost part of the disk, because of the short range of the Lense- 
Thirring effect, while its outer parts tend to remain in its original configuration. The tran- 
sition radius between these two regimes is known as Bardeen-Petterson radius -RbPj shown 
schematically in Figure 1; its exact location depends mainly on the physical properties of 
the accretion disk (Bardeen & Petterson 1975; Kumar & Pringle 1985; Ivanov & lUarianov 
1997; Nelson & Papaloizou 2000; Lubow, Ogilvie & Pringle 2002; Fragile & Anninos 2005). 

The time-scale for alignment can be calculated as (e.g., Natarajan & Armitage 1999): 



where </? is the angle between the black hole spin axis and the direction perpendicular to the 
outer disk. The time derivative of Jbh has the form: 



2. The Bardeen-Petterson eflfect and the alignment time-scale 




(1) 




(2) 
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where it!out is the outer disk radius and L^{r) = Sd(r)nd(r)r^ is its differential angular 
momentum (e.g., Caproni, Mosquera Cuesta & Abraham 2004). is the disk angular 
velocity and is the surface density of the accretion disk integrated over its semi-thickness 
Hfi, derived following Sakimoto & Coroniti (1981) (see also Bardou, Heyvaerts & Duschl 
1998): 

H,{r) = ^^-^^'^^^^^'^ (4) 
where i^nsg = Cs/^d, -f^sg = Cg/(7rGEd) and Cg is the sound speed, defined as: 



^s(r) = A/-r ^^"f-^^^^'^^^^^^'^^^\ (5) 

V a m r a 

where F is the politropic index of the gas, which we have assumed equal to 5/3. 

A rough estimate of Rbp can be obtained comparing the time-scales for Lense-Thirring 
precession and warp transmission through the disk (e.g., Natarajan & Armitage 1999) that, 
on the other hand, will depend on how the warps are being communicated along it. If they 
are transmitted diffusively: 



<P = V^2/^LT, (6) 

where 1/2 is the viscosity acting on the direction perpendicular to the disk; both parameters 
must be calculated at R^p. If the time-scales for the warp and surface density evolution are 
similar 1/2 ~ J^i, where ui is the viscosity along the disk; otherwise 1/2 ~ /(q;)i'i, where f{a) 
is a function of the dimcnsionless viscosity parameter a introduced by Shakura & Sunyaev 
(1973); we adopted f{a) = 2(1 + 7«^)/[«(4 + a^)] as derived by Ogilvie (1999). Assuming 
that the inner radius of the accretion disk is the marginally stable orbit Rms, the viscosity 
can be written as: 



M 



27rEd(r) 



(ilnr2d(r) 
dlnr 



Rn 



(7) 



M is the accretion rate, which is related to the bolometric luminosity Lboi through M / M^dd — 
Lho\/ L-Edd-i where Mgad ^-nd i^Edd ^-re the Eddington mass accretion rate and luminosity, re- 
spectively. 
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In the wave-like regime: 



R^P = c,/Qlt. (8) 



The transition from the diffusive to wave-hke regime occurs at a radius Rt ~ H^/a 
(Papaloizou & Lin 1995). 



3. Application to the accretion disk of NGC 1068 

The physical parameters of the black hole accretion disk system of NGC 1068 were 
determined by Hure (2002) and Lodato & Bertin (2003) from the analysis of the maser hne 
velocities; they are presented in Table 1 and were used in our calculations of the Bardeen- 
Petterson radius and the alignment time-scale. 

The maser velocities present signatures of sub-Keplerian motion; however, in equation 
(3) we used the relativistic Keplerian angular velocity (Abramowicz, Jaroszyhski & Sikora 
1978), since at the radius in which this assumption is not valid, the angular momentum 
becomes negligible. 

We assumed a power- law surface density distribution for the accretion disk E(r) = 
Eo(r/it!g)*, where — GMbu/c^ is the gravitational radius and s = —1.05 (Hure 2002). 
The constant Eq was determined from the mass of the disk M^, derived by Hure (2002) and 
Lodato & Bertin (2003), integrating E(r) from the inner to the outer disk radius. 

The a-parameter was obtained from the expression M — (28.1±0.2)q; Mq yr~^ given by 
Lodato & Bertin (2003). The accretion rate can be calculated from the bolometric luminosity 
Lhoi = ^Mc^, given the accretion efficiency e, which depends on the black hole spin parameter 
a*. We used as lower limit for the bolometric luminosity ~ 7 x 10"^^ erg s~^, as found by 
Gallimore et al. (2004) from the observed free-free emission, and the Eddington luminosity 
-^Edd as upper limit. 

The Bardcen-Petterson radius was obtained from either equation (6) or (8), depending 
if it obeys the diffusive or wave regimes, respectively. 

The calculations were performed for the disk models A and B defined in Table 1, con- 
sidering several values of the black hole spin and the extreme values of the accretion rate 
and a parameter, the latter depending also on the spin value. The results show that the 
Bardcen-Petterson radius for NGC 1068 depends weakly on the accretion rate, being limited 
by the variation of the black hole spin and a parameters between 10~^ and 10~^ pc (about 
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20 and 200 Rg respectively) as shown in left panel of Figure 2. 

For the lower hmit of the a parameter and both black hole masses, Rbp = R^p when 
I a* I > 0.1. For the upper hmit of a, Rbp — R^p for a* < 0.1, turning gradually into R^p 
for larger values of a* at a rate that depends on the a parameter. 

We calculated the alignment time-scale as a function of the black hole spin by solving 
the integral given in equation (3), using as integration limits the Bardeen-Petterson radius 
and an arbitrary outer radius. The results, which are independent of the outer radius value, 
are presented in the right panel of the Figure 2. 

The alignment time-scale turned out to vary between 100 and 10^ yr, while the lifetime 
of the radio jet in NGC 1068 is < 1.5 x 10^ yr (Capetti et al. 1999), similar to the hfetimes of 
AGN activity. Therefore, it indicates that the Bardeen-Petterson effect can perfectly warp 
the inner part of the accretion disk. 

It is important to emphasize that warped accretion disks produced by the Bardeen- 
Petterson effect can be probed by features in relativistically broadened emission iron-line 
profiles (Fragile, Miller & Vandernoot 2005) but unfortunately, since NGC 1068 is a type 2 
Seyfert galaxy, relativistically broadened iron-lines are unlikely to be observed due to the 
obscuring dust torus. 



4. Probing alignment from jet kinematics 

The current angular resolution provided by interferometric techniques is not capable 
to resolve structures with sizes comparable to the Bardeen-Petterson radius in NGC 1068. 
However, the orientation of the inner part of the disk can be traced from the jet, which is 
usually thought (sometimes supported by observations; Ray et al. 1996; Jones et al. 2000) to 
be ejected in the perpendicular direction. In NGC 1068, the jet is not continuous but presents 
discrete features, the upper limit to their expanding velocities being 0.17c (Gallimore, Baum 
& O'Dea 1996). 

The Bardeen-Petterson effect forces the disk to align gradually with the black hole, 
producing precession and a progressive change in the jet direction (Caproni, Mosquera Cuesta 
& Abraham 2004). According to Scheuer & Feiler (1996), the solution of equation (3) gives 
an exponential time variation for the inclination angle between black hole spin and the jet 
direction, as well as for the precession period Pprec^ 



(9) 
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Pprecit) = Poe-^*-*°^/^^"«" (10) 

where (po and Pq are, respectively, the inchnation angle and precession period at time to 
when the disk was formed (to < 0, measured in the past from the present time). Scheuer & 
Feiler (1996) found that the timescales for precession and realignment are identical, implying 
that Pq — Taiign, which will be also consider in our approach. 

We assumed that the jet originated at time to, afterwards each plasma element was 
ejected with constant speed Vjet at a time t, in a direction that forms an angle ip{t) with the 
black hole spin axis, in a plane that rotates with velocity uj{t) = 2tt / Ppi-cc{t) ■ We introduced 
two additional parameters: the viewing and position angle of the black hole spin in relation 
to the line of sight 6 and on the plane of the sky r]. 

In order to compare our model with the observations, we calculated the right ascension 
and dechnation offsets (Act and AS respectively) of each jet element as a function of time 
through: 

Aa{t) ^'^[A{t) cos r] + B{t) sin ri], (11) 

AS{t) = ^ [-A{t) sinr] + B{t) cosr^] , (12) 

where D is the distance to the observer, A{t) = sin ip{t) cos utt and B{t) = sin ip{t) sin urt cos 9+ 
cos ip(t) sin 9. 

We found solutions for several combinations of the input parameters, in the range 25° < 
V^o < 45°, 70° < 6^ < 90° and 14° < < 20°. The parameters to and TaUgn arc scaled by the 
jet velocity ranging between -1.7 x 10^ < to (yr)< -9600 and 7500 < Taiign (yr)< 1.3 x 10^ 
for 0.01c < Vjet < 0.17c (Wilson & Ulvestad 1987; Gallimore, Baum & O'Dea 1996). The 
limits found for TaUgn are perfectly compatible with those derived in last section. 

In Figure 3 we present the comparison between the parscc and kiloparsec scale radio 
maps (Wilson & Ulvestad 1987; Gallimore et al. 2004) and a model with parameters 9 = 80°, 
= 0.17c, to = -9800 yr, Taiign = 7580 yr, ipo = 40° and rj = 17°. We can see that our 
simple kinematic approach reproduces satisfactorily the inverted S-shape of the kiloparsec 
jet, as well as the location of its jet components; at parsec-scales, the position of the jet 
components C and NE, as well as the counterjet knot S2 are also well-reproduced by the 
same model. As it can be seen in the left panel of Figure 3, at larger distances from the 
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core, the full range of position angles provided by our precessing model is not found in 
the observational data, in the sense that the amplitude of the helix towards the northern 
and southern regions are systematically larger than the observed jet aperture. This may 
be due to the implicit assumption of accretion disk rigid-body precession during all the jet- 
time evolution, which might not be totally true, specially close to the initial time to- In fact, 
numerical simulations have shown a period of differential precession preceding the rigid-body 
configuration (Nelson & Papaloizou 2000; Fragile & Anninos 2005). In addition, the galactic 
medium will influence the jet propagation, as can be seen from the [O III] and mid-infrared 
images of NGC 1068 (GaUimore et al. 1996a; Capetti, Axon & Macchetto 1997; GaUiano et 
al. 2003), showing the existence of jet- ambient interaction. 

In all cases considered in this work, the inner disk is not completely aligned with the 
equator of the black hole since </? ~ 11° at the present time. Interestingly, Galliano et al. 
(2003) proposed the existence of a shght tilt of ~ 15° between the Compton thick central 
absorber and the molecular disc in order to reproduce the observed distribution of H2 and 
CO emission in this object. 

5. Conclusions 

We studied the possibihty that the Bardeen-Petterson effect is responsible for the warp- 
ing of the accretion disk in NGC 1068, as derived from the radio continuum and water maser 
observations (e.g., Greenhill & Gwinn 1997; Gallimore et al. 2004). Siich mechanism arises 
from the frame dragging produced by a Kerr black hole whose rotation axis is not parallel 
to that of the accretion disk. 

Using an analytical approach, similar to that suggested by Scheuer & Feiler (1996) and 

Natarajan & Armitage (1999), we calculated the Bardeen-Petterson radius and the alignment 
time-scale between the accretion disk and the equator plane of the black hole for different 
values of the black hole spin. We found that the Bardeen-Petterson radius is roughly limited 
between 10~^ and 10"'^ pc, while the alignment time-scale ranges from about 100 to 10^ yr. 
Those estimates are perfectly compatible with the upper limit to the AGN lifetime, derived 
from the jet kinematics. 

We also showed that the general form of the parsec and kiloparsec scale jets can be 
reproduced by the model in which the jet is precessing with a period Pprec and aperture 
angle ip that decrease exponentially with a time-scale Taiign, as expected from the Bardeen- 
Petterson effect. 
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Fig. 1. — Schematic view of the Bardeen-Petterson effect. An accretion disk with inner and 
outer radius Ra,s and i?out respectively, having a semi-thickness and misahgned initially 
by an angle (/? in relation to the angular momentum of the black hole Jbh and mass Mbh, will 
be warped by the Bardeen-Pcttcrson effect. The Bardccn-Petterson radius i?BP marks the 
transition between the aligned and misaligned disk in relation to the black hole's equator. 
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Fig. 2. — Left panel: Bardeen-Petterson radius as a function of the black hole spin for 
NGC 1068. Each point corresponds to the mean value of Rbp considering the limits of 
minimum and maximum accretion rate. Circles and triangles refer respectively to the models 
A and B. Filled symbols are related to the lower value of a-parameter, while open ones to 
its upper value. Right panel: Time-scale for alignment between accretion disk and the black 
hole's equator, using the same nomenclature for the symbols. 
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Fig. 3. — Snapshot of the changes in the jet orientation due to the Bardeen-Petterson effect. 
Left panel: Kiloparsec scale jet of NGC 1068 at 5 GHz (Wilson & Ulvestad 1987). Solid hne 
indicates the model prediction from equations (9) and (10), assuming 9 — 80°, Vjet — 0.17c, 
to — —9800 yr, TaUgn = 7580 yr, cpo — 40° and rj — 17°. Right panel: Comparison between 
the map at 1.4 GHz (Gallimore et al. 2004) with the same model shown in the left panel. 



Table 1. Physical parameters of the disk. 



Parameter 


Model A 


Model B 


Mbh (Mq) 


(1.2 ± 0.1)xlO'' 


(8.0 ± 0.3)xl0« 


s 


-1.05 ± 0.10 


-1.05 ± 0.10 


Md(r=lpc) (Mq) 


(9.4 ± 1.6) xlO^ 


(8.6 ± 0.6) xlO** 


So (g cm-2) 


(1.06 ± 0.23) X 10'' 


(1.49 ± 0.18) xlO'^ 


LEdd (erg s-i) 


(1.5 ± 0.1)xl0-*5 


(1.0 it 0.3) x 10*5 


M/Mjsdd 


0.46 - 1.0 


0.71 - 1.0 


ae(a*) 


(4.4- 9.4) xlQ-" 


(4.4 - 6.3) X 10"* 



Note. — The quoted errors in Mbh and s are given by Hure 
(2002) and Lodato & Bertin (2003), while for Sq were obtained 
from error propagation. 



